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We consider theoretically a superconducting qubit - nanomechanical resonator (NR) system, which was re- 
alized by LaHaye et al. [Nature 459, 960 (2009)]. First, we study the problem where the state of the strongly 
driven qubit is probed through the frequency shift of the low-frequency NR. In the case where the coupling 
is capacitive, the measured quantity can be related to the so-called quantum capacitance. Our theoretical re- 
sults agree with the experimentally observed result that, under resonant driving, the frequency shift repeatedly 
changes sign. We then formulate and solve the inverse Landau-Zener-Stiickelberg problem, where we assume 
the driven qubit's state to be known (i.e. measured by some other device) and aim to find the parameters of the 
qubit' s Hamiltonian. In particular, for our system the qubit's bias is defined by the NR's displacement. This 
may provide a tool for monitoring of the NR's position. 



I. INTRODUCTION 

Nanoelectromechanical systems have recently attracted at- 
tention because of both possible applications (e.g. in sensing) 
and interest in fundamental quantum phenomena in meso- 
scopic systems.^ Particularly interesting is the coupling of the 
mechanical motion of a nanomechanical resonator (NR) to an 
electric mesoscopic system. A few examples are carbon nan- 
otube NRs coupled to electron transport^ and a metallic NR 
coupled to an LC tank circuit?. It was proposed theoreti- 
cally that for sensing and controlling the NRs, superconduct- 
ing few-level circuits (qubits)"^ can be effectively used.^'^ For 
example this approach was applied in the demonstration of 
the ground state of a high-frequency piezoelectric dilatational 
resonator coupled to a superconducting phase qubit. ^ 

Successful coupling of a NR (a suspended silicon nitride 
beam) to a charge qubit allowed LaHaye et al. |8] to demon- 
strate both ground- state measurement and excited- state spec- 
troscopy as well as Landau-Zener-Stiickelberg (LZS) interfer- 
ometry of the qubit. The spectroscopy was performed with 
weak driving, where the position of the resonance gave the 
information about the qubit levels. In the regime of strong 
driving, where the qubit's evolution experiences repeated LZS 
transitions at the avoided crossing, the resulting interference 
is visualized in the LZS interferograms |9]. The LZS interfer- 
ometry was demonstrated on superconducting qubits probed 
by different methods (see Ref. [9] and references therein), 
as well as studied for other different physical realizations of 
strongly-driven two-level systems in Refs. QO]. 

In the work by LaHaye et al, Ref. JH], the NR's frequency 
shift was used for monitoring the qubit's state. For the theo- 
retical description of the NR-qubit system, the perturbation- 
theory procedure developed in Ref. |I5|] was used. The theory 
says that the NR's frequency shift Acjnr is negative for the 
qubit in the ground state and zero when the two qubit states 
are on average equally populated under the periodic driving. 
This allowed to describe the ground-state and low-amplitude 
spectroscopy measurements.^ However, this theory does not 
explain the experimentally observed sign changes of Acjnr in 
the strong-driving regime, where the frequency shift becomes 



positive. 

In this work we consider the NR-qubit system semi- 
classically. Within this approach, we describe the qubit as 
a quantum system coupled to a classical resonator, with the 
oscillation-energy quantum much smaller than the thermal en- 
ergy, fic^NR <C ksT. Note that such a semi-classical approach 
was successful for the description of most phenomena related 
to atom-light interaction.^^ 

The impact of the qubit on the resonator's frequency shift 
can be described in terms of the so-called quantum capaci- 
tance, as studied for the qubits in Refs. [12,131. The quantum 
capacitance is defined as the derivative of the average charge 
on the qubit with respect to the applied voltage. The charge 
can then be related to the charge-qubit occupation, the deriva- 
tive of which (under resonant driving) exhibits sign changes. 
Similar sign-changing response under strong driving was re- 
cently studied for qubits probed by an LC (tank) circuit for 
capacitive coupling ^"^-^^ as well as for inductive coupling 
Thus, in the first part of this work (Section II) we study the 
situation where the strong-driving qubit's state is probed by 
the NR. 

In Section III, we formulate the inverse problem. There, we 
are interested in the influence of the NR's state (its position) 
on the qubit's state. We graphically demonstrate the formula- 
tion of the problem for the direct and inverse interferometry 
in Fig.[T] There, the two-level system represents a qubit with 
control parameter sq, the parabola represents the resonator's 
potential energy as a function of the displacement x. Thus, in 
the first part of our work (Sec. II) we deal with the direct prob- 
lem, where the influence of the qubit's state on the resonator 
is studied. 

The second part of this work (Sees. Ill and IV) is devoted to 
the inverse problem, where we study the influence of the res- 
onator's state on the qubit's state. Measuring the latter is an 
alternative method for defining the NR's displacement. This 
approach can be related also to other inverse problems for 
two-level systems, as studied in Refs. |[T8| - [2Q|] . Generalization 
of the results can also be applied to other quantum systems for 
which the problem of defining the Hamiltonian 's parameters 
with given system's state was studied in Ref. ll2lll . In Section 
IV we demonstrate how the inverse problem can be solved for 
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FIG. 1: (Color online) Schematic representation of the formulated 
problems for direct and inverse interferometry. The red curves on 
the left represent the bias-dependent energy levels of the qubit, and 
the green parabola on the right shows the potential energy of the 
(classical) resonator. In the direct problem, the resonator is used to 
probe the state of the qubit. In the inverse problem, the response of 
the qubit to external driving is used to infer the state of the resonator. 



different driving regimes in a generic two-level system, and 
we comment on the possibility of applying this technique for 
superconducting qubit-NR systems. 



II. CHARGE QUBIT PROBED THROUGH THE 
QUANTUM CAPACITANCE 

The split-junction charge qubit (also called Cooper-pair box 
and shown in red in Fig. [2]) consists of a small island between 
two Josephson junctions. The state of the qubit is controlled 
by the magnetic flux <l> and the gate voltage Vqpb + ^mw- 
Here Vcpb is the dc voltage used to tune the energy levels 
of the qubit and Vuw = sin ujt is the microwave signal 
used to drive and manipulate the energy-level occupations. 
The Cooper-pair box is described in the two-level approxi- 
mation by the Hamiltonian in the charge representation (see 
e.g. Ref. 1 8] and Appendix A) 



^0 



A sin ujt 



(1) 



Here the tunnel splitting A is equal to the Josephson energy 
^j, which is controlled by the magnetic flux ^ 

A = ^j = ^jo|cos(^^/^o)|. (2) 
The charging energy and the driving amplitude are given by 

£0 = 8^c(ng - 1/2), A = SEcn^, (3) 

where the Coulomb energy Eq = 6^/2(7$] is defined by the 
total island's capacitance Cy: = 2Cj + Ccpb + C^r, de- 
fined with the notation 2Cj = Cji + Cj2; the dimension- 
less driving amplitude is = CcPBVi^/2e; the dimension- 
less polarization charge Ug = tinr + ^cpb is the fractional 
part of the respective polarization charges in the plates of 
two capacitors: tinr = {A/'nr} and ncpB = {^cpb} with 
N^R = C^RV^R/2e and A^cpb = C'cPBl^cPB/2e. 

Here we consider the Cooper-pair box formed by four ca- 
pacitances, Cji, Cj2, CcpB, and Cnr (Cj > Ccpb,C'nr). 
One of the plates of the latter capacitor is formed by the NR, 
which is characterized by the displacement at the midpoint x. 
This displacement is usually much smaller than the distance d 
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FIG. 2: (Color online) Schematic diagram of a split-junction charge 
qubit coupled to a nanomechanical resonator. The charge qubit 
(shown in red) is biased by the magnetic flux $ and the dc+fiw volt- 
age, Vcpb + Vmw, to which it is coupled through the capacitance 
CcPB. The qubit is coupled to the NR (shown in green) through the 
capacitance Cnr. The NR is biased by a large dc voltage Vnr; its 
state is controlled and measured by applying the dc and rf voltages 
between the gate and the NR, Vqnr and Vrf, through the capaci- 
tance Cgnr. The NR's motion is described by the displacement at 
the midpoint x. Capacitances form the island (Cooper-pair box) with 
the total capacitance Cs, voltage Vi and charge —2en. 



between the plates, in which case the capacitance between the 
NR and the qubit reads^^ 



C'nr(^) ^ C'nro + 



NR 



dx 



Cnro 1 



1 dQ 



NR 



Cnro dx 



(4) 



(5) 



(By the subscript here we mean the values at x = 0; in 
what follows this subscript is assumed). The displacement of 
the NR influences the qubit through the changes in the polar- 
ization charge; to make this influence significant, a large dc 
voltage Fnr (of the order of volts) is applied. On the other 
side, the NR is biased by dc and rf voltages, Vqnr and Vrf, 
through the capacitance Cqnr, which provide its control and 
read-out 

The influence of the qubit 's dynamics on the nanomechan- 
ical resonator can be described in different ways. In Ap- 
pendix A we present a detailed derivation of the influence of 
the qubit 's state through the voltage Vi and the average polar- 
ization charge — 2e (n) of the CPB on the NR's dynamics. An 
alternative, and maybe physically more illustrative, approach 
is to describe the CPB as an effective capacitor, which is the 
subject of Appendix B. Here, in the main text, we present only 
essential results, referring the interested reader to the Appen- 
dices. 

As a result of the interaction between the qubit and the 
NR, the resonance frequency of the NR is shifted (see Ap- 
pendix A). The result can be written in the following form 
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The frequency shift Acjnr is defined by the derivative of the 
average extra Cooper-pair number on the island (n) = • Pq + 
1 • Pi = Pi. Here Pq (Pi) stands for the probabiHty of having 
0(1) extra Cooper pair. 

Alternatively to the approach above, the effect of the qubit 
on the NR can be described in terms of the effective (dif- 
ferential) capacitance, as described in Appendix B, Ceff = 



dQ^^/dV, 



NR 



Cg- 



Cq, where the relevant quantum 



capacitance is given by 



Cc 



drig- 



(8) 



The term "quantum" capacitance is used here to denote the 
(small) qubit- state-dependent addition to the classical (ge- 
ometric) capacitance. Obviously, Eq. (|6]) can be rewritten 
in terms of the quantum capacitance (cf. discussion in Ap- 
pendix C for the qubit-LCP circuit system) 



NR 



^NR 



(9) 



of the resonance, to negative, to the right of the resonance 
point. Thus, the resulting behavior of the observable (either 
Acjnr or Cq) is defined by the competition of the two terms 
in Eq. (fTTl) . In what follows we will use Eq. (fTTl) for the su- 
perposition states (which appear under driving). Note that 
a similar approach for calculating the effective (quantum) in- 
ductance was used in Refs.[T6lfT7l 

The dissipative dynamics can be described with the Bloch 
equations written in the energy representation (where relax- 
ation appears naturally). To characterize dissipation we use a 
result of the spin-boson model with the spectral density de- 
fined with the dimensionless parameter a, J{uj) = ahuj, see 
e.g. Ref. ^25 and references therein, while the low-frequency 
1// noise is described by the peak of J{uj) at ^ 0. Then 
the relaxation and dephasing times are defined by the spectral 
density at a; ^ AE and u ^ respectively as following 
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where /3 = (Cs/Cnr)/3. 

The qubit' s density matrix in the energy representation (in 
the eigenbasis of the time-independent Hamiltonian) can be 
parameterized in terms of the respective Pauli matrices r^, p = 
^ (Xr^ + Yry ^Zr^), as e.g. in Ref. (TtI]. Here Z = (r^) 
is the difference between the occupation probabilities of the 
excited and ground states. Now we express the probability 
of having one excess Cooper pair. Pi , by changing from the 
energy basis to the charge basis, and obtain 



Pi = ifl- Ax + 4^Z), A^ 
^ 2 V AE AE ^ ' 



A^^el (10) 



And this gives (after time- averaging over the driving period 
27r / uj) for the quantum capacitance the following 



Cq 



dZ 



AE^ 



2AE dn. 



(11) 



where we have taken into account that in the stationary state 
X averages to 0.^ 

As we can see from Eq. (fTTt . the quantum capacitance is 
defined by the value Z = (tz). In particular, we obtain the 
quantum capacitance and the respective frequency shift in the 
ground/excited (g/e) state with Z = ±1 
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(12) 



This result, obtained in the semi-classical approach, is in 
agreement with the one obtained in Ref. Is] and used in Ref. [H. 
Equation (fTTl) is a more general result, where the second 
term describes the sign-changing behavior near resonance. 
Namely, when sweeping the gate voltage ng, the quantity Z 
changes from —1, far from resonance (in the ground state), to 
in resonance (when the levels are equally populated). This 
describes the maximum of Z in resonance and the change 
of its derivative dZ/drig from positive, in the left vicinity 
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Here the (relatively large) phenomenological parameter B 
was introduced to describe the low-frequency 1// noise. We 
note that alternatively the low-frequency noise could be taken 
into account as the averaging of the final solution resulting in 
some blurring of the resonances, as e.g. in Ref. 14. The values 
for the relaxation and dephasing times define the shape of the 
resonances (as for example it is later described by Eqs. ([28l 
[31])). In this way, the width of the resonances can be used 
for the estimation of the dephasing rate. In our case, we have 
taken a and B as the fitting parameters, to obtain better re- 
semblance with the experimental results. 

We display the direct LZS interferometry in Fig. [3l where 
the resonator's frequency shift Acjnr was calculated with 
Eqs. ([9]) and (fTTl) . Figure [3] demonstrates that our formal- 
ism is valid for a description of the experimentally measurable 
quantities: the quantum capacitance or the resonant frequency 
shift^'^"^ (see also Appendix C). Such a description allows 
to correctly find the position of the resonance peaks in the 
interferogram and to demonstrate the sign-changing behav- 
ior of the quantum capacitance, which relates to the measur- 
able quantities. The appearance of the interferogram depends 
on several factors: the values of the qubit parameters, the 
model for the dissipative environment (such as Eqs. ([T3l fT4l) 
and the parameters a and P), the value of the bias current 
(which distorts the shape of the resonances, as demonstrated 
in Ref. |17]). Moreover, the formalism presented above is 
valid for the case where the qubit 's dynamics is much faster 
than the NR's dynamics; otherwise one should study the co- 
operative dynamics of the composite system; see, e.g., discus- 
sions in Refs. 1 14] and 1 17]. However, we will not go here into 
more detailed calculations, since our aim was to demonstrate 
the simplest approach for the description of the experiment in 
Ref. 1 8]. 
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FIG. 3: (Color online) LZS interferometry probed via the resonator's 
frequency shift Acjnr. (a) The frequency shift versus the energy 
bias (rig) and the driving amplitude (ri/x). Arrows show the values 
of ri/j^ and rig at which the graphs (b) and (c) are plotted as func- 
tions of Ug and n^, respectively. The upper curves were shifted for 
clarity. The parameters for calculations were taken close to the ones 
of Ref. 1 8]: cJNR/27r = 58 MHz, Ejo/h = 13 GHz, Ec/h = 14 
GHz, cj/27r = 4 GHz, ksT/h = 2 GHz, a = 0.005, B = 0.2, and 
the proportionality coefficient [3 defined by the qubit-NR coupling 
constant A from Ref. [8]: hX^ /ttEjq = P • Ec(jO]sir /ttE jo = 1.6 
kHz. 



ways to probe the qubit's state). Given the known qubit state, 
we aim to find the Hamiltonian's parameters. Particularly, we 
are interested in the parameter-dependent bias eo{x). 

On one hand, we can study here the general ("reverse engi- 
neering") problem in the spirit of Refs. [18.191. On the other 
hand, we aim to provide the basis for measuring the NR's po- 
sition X by means of probing the qubit's state, while x = x{t) 
is considered a slow time-dependent function. 

In what follows we will consider the driven qubit's state 
with emphasis on finding optimal driving and controlled offset 
parameters (A, uj, and Eq) for the resolution of the small bias 
component Ssq . We will assume that the dynamics of the pa- 
rameter X is slow enough not to be considered during the mea- 
surement process. Depending on this slowness, the measure- 
ment might have to involve only one passage of the avoided 
crossing, or it can involve long-time driving and stationary- 
state equilibrium of the qubit. Our aim is to find a sensitive 
probe for small Seq. For high sensitivity we require substan- 
tial changes in the qubit's state for small changes of Eq given 
by Sso- For a quantitative definition of the sensitivity one can 
consider the derivative of the probability with respect to the 
bias £o. 



IV. RESULTS FOR THE INVERSE LZS 
INTERFEROMETRY 

In this section we consider the inverse problem for the 
qubit's dynamics, in particular how to infer the qubit's bias £o 
from the measured qubit state. For concreteness, we consider 
the qubit driven by the bias e{t) = eo-\-A sin ut. For purposes 
of analyzing the short-time dynamics, one would consider a 
single passage or a sequence of a small number of passages 
through the avoided level crossing. If the time-dependence of 
the bias Eq {x) is so slow that the multiple-passage dynamics 
is relevant, then the stationary qubit state can be considered. 



III. THE BIAS INFLUENCED BY THE RESONATOR: 
PROBLEM FOR THE INVERSE INTERFEROMETRY 



A. Single passage: non-linearity in the Landau-Zener problem 



Let us now consider the qubit's bias Sq, Eq. ©, as a func- 
tion of the NR's displacement x. For small x <C ^, we have 
the expansion d?]), which results in the decomposition of the 
bias 



where 



eSK) = 8Ec(ng-l/2), 

X 

Seo{x) = S^^c^NR^. 



(15) 

(16) 
(17) 



Here we have used the fact that x <C ^ and Cnr ^ C^] • 

The Hamiltonian of the qubit ([T]) with the parameter- 
dependent bias eo{x) brings us to the following problem. Let 
us assume that the qubit's state is known (i.e., this is measured 
by a device whose details we do not consider here for simplic- 
ity; see Refs. |[T2|fT3lfT6l|230 for different realizations of the 



The linearization of the bias in the vicinity of the 
avoided crossing (where e{t) = 0) results in the approx- 
imation that this region is swept at the -dependent rate 
Auy^l — (eo/A)'^ (for details see Ref. [9]). The respective 
probability of the non-adiabatic transition to the upper adia- 
batic level is given by the Landau-Zener formula 



Plz exp 



2 Ahu 



(18) 

In other words, the non-linear dependence of the bias on time 
has the effect that the Landau-Zener probability depends on eq 
(see also Ref. H^]), which is demonstrated in Fig.^a). We 
note that here | | < A and the formula ([TSl) gives numerically 
incorrect results when sq tends to A. 

To quantify the sensitivity of the transition probability to 
small changes in the bias, in Fig.|4tc) we plot the derivative of 
the excitation probability P^"^ with respect to Sq- We can see 
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FIG. 4: (Color online) Upper-level excitation probability P+ after 
(a) single passage and (b) double passage, plotted for A/ A = 5 and 
/kj/A = 0.2, versus the bias so- The sensitivity to the changes of 
the bias so, defined as the derivative, x — \dP+ / d{eo / A)\, is plotted 
in (c) and (d), respectively. Solid lines were plotted with Eqs. ([Tsl 
and (I2TI) . while dashed lines were calculated numerically. 



anticrossings at t2 and ti + 27r/a;: 

C2 = ^y^ ^/^TTWdt, (22) 

and (fs is the Stokes phase. 

Stiickelberg oscillations, described by Eq. ([2T]) . are demon- 
strated in Fig. lUb) for < Eq/A < 1. The respec- 
tive sensitivity is shown in Fig. HJd). The agreement of 
the analytical formulas and numerical calculations is remark- 
able (as demonstrated in Fig. |4|. One can notice that the 
sharper the Stiickelberg oscillations, the higher the sensitiv- 
ity. This is related to the period of the Stiickelberg oscilla- 
tions, which decreases with increasing A/cu. Here we also 

note that P\^^\so) is not a symmetric function, and the pe- 
riod of the Stiickelberg oscillations is smaller for Sq < than 
for £0 > 0. Therefore, using negative values of Eq results in 
slightly higher sensitivity than what is shown in Fig.|4td). 

The factor Plz(1 — Plz) in Eq. (|2T]) is described by the 
one-passage problem above. Consider the term cos^ ^2- For 
£q = and Seq <C A we have^ C2 ^ ~ f example, 
for ^ = 2/c7r + ^ we obtain 



that the non-linearity of the bias results in an increase of the 
sensitivity. 

For the single-passage case it is straightforward, from 
Eq. (fTSl) . to find the solution for the inverse problem eq = 
eo{P\^'^)- In particular, in the case Sq = and 5eo <C A we 
have 



LZ 



LZ,0 



7 f feo 
2 \ A 



Plz,o 



(19) 



and the solution for the inverse problem becomes 




(20) 



B. Double passage: Stiickelberg oscillations 

Next, consider the situation where the avoided crossing re- 
gion is passed twice. For example, the qubit can be driven 
by a sinusoidal pulse of length 27t/uj. Alternatively, triangu- 
lar pulses can be used to drive the qubit twice through the 
avoided-level crossing, as in Refs. iHM. In both cases, the 
double-passage process can make use of quantum interference 
to increase the sensitivity of our problem through the accumu- 
lation of the Stiickelberg phase. 

The upper-level excitation probability after the double- 
passage is^ 



f^)^2PLz(l-PLz)(l + ^^ 



(23) 



This describes a linear dependence on the small bias Seq, 
which is a significant increase in sensitivity as compared to 
the quadratic dependence on Seq in the single-passage case 
above, Eq. ([T9l) . If the decoherence is negligibly small, one 
can further increase the sensitivity of the excitation probabil- 
ity to small changes in the bias due to interference by consid- 
ering multiple-passage case. 

The formula ([23l) can be conveniently used to make quan- 
titative estimates. Consider this for the example of the qubit- 
nanomechanical resonator system as in Ref. | 8]. First, to in- 
crease the sensitivity of the changes of P^^"^ with respect to 
Sso, we choose the smallest possible frequency u. In our 
case the driving period should exceed the decoherence time 
T2 and the NR oscillation period 27r/co'NR- For superconduct- 
ing qubits T2 is typically higher than 1 /is. Then, we are 
limited by the relation uj > cjnr, and we take uj/27t ~ 0.1 
GHz. We choose the parameters A{n^) and A(<l>) such that 
Plz ^ 1/2. Assuming nNR = I and 8Ec/h = 100 GHz, we 
obtain the change of the probability with changes in the NR's 
displacement AP^^"^ = lO^x/^. This means that for probing 
a displacement ofx ^ 10~^^, one has to be able to measure 
population changes p!^^-' ^ 0.01. This level of accuracy is 
achievable with superconducting qubits.— 



C. Multiple passage: stationary solution 



Pi''^ = 4Plz(1 - Plz) sin'(C2 + ^s). 



(21) 



where (2 is the phase acquired during the evolution between 



Now we assume that what is relevant for our inverse prob- 
lem is the stationary state of the driven qubit. To analyze the 
analytical expressions, we consider two limiting cases. 



6 




FIG. 5: (Color online) Slow-passage and fast-passage LZS interferometry of a qubit. (a) and (d): the time- averaged upper-level occupation 
probabilities, defined in the adiabatic (P+) and diabatic (Pup) bases, as functions of the bias eo and driving amplitude A. The parameters are 
the same as for Fig. [3] except for the frequency: (a) uj^i^ — 6.5 GHz < /\/h and (d) uj/2n = 20 GHz > A/h. (b) and (e): Cross-sections for 
the respective dependencies of the upper-level occupation probabilities as functions of the bias along the horizontal dashes shown in red and 
green in (a) and (d). (c) and (f): Inverse graphs, which show the dependence of the bias on the upper-level occupation probabilities (assuming 
that So lies on the right-hand side of the resonance peak). 



1. Slow -pas sage limit 



we have 



For the analytical description of the upper-level occupation 
probability in the adiabatic limit, when 7 > 1, we use the 
following formula from Ref . t90 



where 



Plz(1 -cosC+cosC-) 
^ ~ sin^ C; + 2Plz (1 - cos cos (- ) ' 

c; = C1 + C2, c- = Ci-C2, 



(24) 



(25) 



(26) 



Analyzing the interferogram in Fig. [3a), we find the possi- 
bility to obtain a sensitive working point with a driving am- 
plitude a little bit lower than the one where the width of the 
resonance line tends to zero, that is 2A/huj = 2iTn — a, a <C 1 
[see the red and green dashes in Fig.Oa)]. It follows that 



-j^ 1 ^LZ {TrSeo/huj) 



2a2 + PLz (^feo/M 



2 ' 



(27) 



and (2 is given by Eq. (1221) . Formula (l24l) is illustrated in 
Fig- Ha). Consider = 0, then for strong driving, A A, 



which is equal to zero at Seq = and quickly tends to 1 /2 
with increasing dso. This is demonstrated in Fig.[5tb). 



2. Fast-passage limit 

In the fast-passage and strong-driving regime (where 7 <C 
1), the rotating- wave approximation gives for the upper-level 



occupation probability^^ 



k T1T2 

A Jfe (A/fuj) 



(28) 



(29) 
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where Jk is the Bessel function. Formula (|281) is demonstrated 
in Fig. [3d). If the relaxation is not taken into account, then 
in the vicinity of the k-\h resonance (where Eq = khu) we 
obtain the Lorentzian dependence on the small bias shift Seq : 



^ 1 



up 



2Ssl 



(30) 



This describes the resonance peak, P^p = 1/2, at Ssq = 0, 
which is demonstrated graphically in Fig. Oe). Its width is 
defined by and is minimized for the values of A/hu in the 
vicinity of the zeros of the Bessel function. With relaxation 
taken into account, the sensitivity is defined by the half- width 
of the resonances, given by 



T2 



(31) 



This means that to increase the sensitivity, which is related 
to the sharpness of the resonances, one has to decrease the 
decoherence rate. 

Here we note that it was assumed that the measurement 
time is much smaller than the resonator's period, T^eas <^ 
27r/cjNR- On the other hand, to reach a stationary state, the 
measurement time should be larger than the relaxation time, 
Pi, 2 < Pmeas- This mcans that the results presented in this 
subsection are relevant for qubits with short relaxation times 
and for resonators with small frequencies. Alternatively, one 
should solve the problem which explicitly takes into account 
X — x{t). 

Formula ([3T1) allows us to make estimates, as we did at the 
end of the previous subsection. For A/huj equal to one of 
the Bessel-function zeros and for T2 = 4ns <C 27r/a;NR, we 
obtain that the probability Pup changes by about 1/4 when 
the bias changes by Aeo/h ~ 0.25 GHz. On the other hand, 
we have seen that Sso/h ^ 100(x/^) GHz. This means that 
in order to observe changes x ~ 10~^^, one has to distinguish 
changes in P^p ~ 10 which is also possible, in principle. 



D. Inverse interferometry: qubit probes resonator 

The idea of the measurement procedure, presented in Fig.[5l 
could be as follows. Driving the qubit in a wide range of pa- 
rameters is done first to plot the interferogram as in Fig. Oa) 
and/or (d). Then a region of high sensitivity, where small 
changes in the qubit bias result in large changes in the final 
state, is chosen. Examples of such high- sensitivity regions are 
shown in Fig. Ob) and/or (e). 

From Fig. [5] we can see that both the slow-passage limit, 
demonstrated in Fig. Oa-c), and the fast passage limit 
[Fig - Hd-f)] can be used for the solution of the inverse prob- 
lem. The choice of the optimal working point and its vicinity 
will depend on the specific parameters of the problem. For 
illustration, in Fig. Oa) and (d) we marked by red and green 
small dashes two possibilities of having the dip in Fig. Ob) or 
the peak in Fig. Oe) being narrow (red curves) or relatively 
wide (green curves). 



In principle, a low-amplitude slice near the bottom of 
Fig. Od) can be used to obtain a sharp resonance peak, as 
in Fig. Oe). However, based on the results of Refs. 11911301] . 
it seems that the width of the resonances might be increased 
more for low- amplitude driving due to the influence of the 
noise and decoherence. From the experimental point of view 
the best strategy is probably to obtain a wide range interfero- 
gram and then choose a narrow resonance. 

One can now bias the qubit at a high-sensitivity point, apply 
a "measurement pulse" to the qubit, measure its state at the 
end of the pulse and extract the resonator's position x from the 
measured qubit's state, see Fig. [3c) and (f), where eo (which 
parametrically depends on x) is plotted as a function of the 
qubit's occupation probability. 

It should be noted here that the measurement pulse, which 
is essentially a driving signal applied to the qubit, can take a 
short duration at the beginning of the measurement process. 
Afterwards the final state of the qubit is read out in the ab- 
sence of any driving fields. As a result, issues that only affect 
the qubit on relatively large timescales, e.g. dephasing and the 
slow measurement of the qubit's state, do not affect the qubit's 
ability to accurately measure the instantaneous position of the 
resonator. It should also be noted that this measurement pro- 
cedure is a single- shot type of measurement and not a contin- 
uous measurement. One could in principle use several qubits 
in order to perform multiple measurements on the state of the 
resonator. 



V. CONCLUSIONS 

We have analyzed a measurement scheme where a qubit 
is probed via a quantum capacitance. We demonstrated the 
sign-changing behavior of the quantum capacitance where the 
strongly-driven qubit exhibits a LZS interferogram. Our semi- 
classical calculations were used to describe recent experimen- 
tal results^ for the LZS interferometry of the qubit probed by 
a NR. 

Then, motivated by the experimental work by LaHaye et 
al. 1 8], we formulated the inverse problem. The inverse LZS 
problem was formulated and solved for a generic two-level 
system in several driving regimes. More specifically, we have 
split the quasi-constant bias Eq into an externally-controlled 
part SQ{ng) and a small part S£o{x) that is to be measured 
through the qubit's state. For the qubit-NR system the former 
can be changed through the gate voltage to realize the most 
efficient measurement working point; the latter was assumed 
to be a function of the NR's displacement x. 

We have shown how the inverse problem can be used for 
defining the NR's displacement. First, one should find (mea- 
sure) the direct LZS interferogram (in a wide range of parame- 
ters). This allows finding the qubit's parameters and choosing 
the optimal bias Eq. Then, fixing the qubit's parameters at the 
optimal working point, small changes due to the slow NR's 
motion may be used for measuring its displacement. 
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Here we note that to obtain the charging Hamiltonian of 
the CPB in the two-state approximation, we consider A^g = 
N -\-ng close to a half-integer number, where N is the integer 
part of Ng, and Ug = {Ng} is the fractional part. Then, with 
n = N -\- n and < 1, we obtain for i^cPB = 
the charging part of Hamiltonian, Eq. ([T]). Here the operator 
for the extra Cooper-pair number n = {1 -\- az)/2 acts on the 
"charge" basis states as follows: n |0) = and n |1) = 

At this point we assume that the qubit's dynamics is much 
faster than that of the classical NR, so the equation for the 
NR can be averaged over the period 27t/uj and then the NR's 
dynamics is defined by the time-averaged voltage 



Appendix A: Semi-classical theory for the qubit-resonator 
system 

In this Appendix we consider the semi-classical theory for 
the qubit-NR system. The equation for the displacement x 
of the classical NR with effective mass m, quality factor Q, 
eigenfrequency ujq and driven by the external force F, is 



m 



X muo dx 



Q dt 



mujnX 



F. 



(Al) 



In our problem, presented in Fig. O the NR is influenced by 
the voltage difference from both sides. On one side (to the 
right of the NR in Fig. O the voltage difference contains the 
large constant part, AV = V^r — Vqnr, and the small rf 
driving component, Vrf = Va cos uj^ft. The force due to 
these voltages is 



^GNR 



1 fdCcNR] 

2[ dx J 



Vgnr — ^rf)^] 

-FACOSCJrf^, (A2) 



where Fa = {dCcNR/dx) • AV • Va- From the other side 
(left side of the NR in Fig.O the voltage difference is defined 
by the island's voltage Vi. The respective force is 



^NR 



^^[Cnr(Fnr 



NR 



dx 



"^NR^ (C'nR^i) 



(A3) 



In the Coulomb-blockade regime, the voltage Vi is defined 
by the quantum-mechanically averaged island's charge — 2en, 
which is given by the sum of the charges on the plates of the 
capacitors that define the island, 

- 2en = Qji + Qj2 - Qcpb - Qnr- (A4) 
For the island's voltage it follows that 



Vi 



2e{Ng + sinujt — n) 



(A5) 



C'nr^nr , C'cpbVcpb 



2e 



2e 



A^NR + iVcPB. (A6) 



- 2e(ng - (n)) 
"^^^^ • 



(A7) 



In what follows this time averaging is assumed. 

Denoting the sum of the constant terms in Eqs. (IA2IIA3I) as 
Fo, we obtain 



dF 

F = Fq ^—x - Fa coscjrft, 
ox 



OF _ 2 / CnrFnr V 



d{n) 



(A8) 



(A9) 



The term Fq results in an (irrelevant) constant displacement 
of the NR, while the linear term results in the resonance fre- 
quency shift in Eq. (lAlb as follows 



2 dF ^2 

Then we obtain the NR's frequency shift 

1 dF 



(AlO) 



Auj- 



NR 



• ^NR — ^0 



2muJo dx 



Auji^Au2, (All) 



where Aoui and Acj2 correspond to the two terms in Eq. 
The term Aui does not depend on the state of the qubit; we 
therefore define the qubit- state-dependent frequency shift 

A'^NR = A^NR — A<^1 = (A12) 

which leads to Eq. Q. 

Appendix B: Quantum capacitance 

In addition to the theory presented in the previous Ap- 
pendix, it is useful to consider the system qubit-resonator by 
introducing the quantum capacitance, which is the subject of 
this Appendix. 

Let us introduce the effective (differential) capacitance, as 
it is shown in Fig. Oa), by differentiating the charge Qnr 
of the capacitance Cnr as follows^^: Ceff = 9Qnr/S^nr- 
Then, for the charge Qnr = (^nr — V'i)Cnr with the island's 
voltage given by Eq. (IA7I) , we obtain 



C'eff — Cg 



(Bl) 
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i C 



CPB 



i 2C, C 



NR 



NR 



(b) 





(c) 



FIG. 6: (Color online) Scheme showing how the charge qubit can be 
described as an effective capacitance coupled either to the NR or to 
LCR resonator, (a) To the left, the charge qubit (CPB) is shown to 
be described as the capacitance 2Cj controlled by the voltage Vcpb 
and coupled through the coupling capacitance Cnr to a measuring 
circuitry. This is described as the effective capacitance Ces as shown 
to the right, (b) The effective capacitance is coupled to the NR, which 
can be used to model our system shown in Fig. [2l (c) The effective 
capacitance is coupled to the electric LCR tank circuit. 



which consists of the quantum capacitance Cq, given by 
Eq. (O, and the geometric capacitance Cgeom 



Co-, 



2Cj + Cnr 



(B2) 



where the latter approximation is valid for Ccpb ^ Cj , Cnr. 

Alternatively to the approach of the previous Appendix, one 
can consider the force Fnr as the electrostatic force from the 
effective capacitance [see Fig. Ob)]: Fnr = {^ei^V^^). 
Then the term with the quantum capacitance, in which ^ 

^NRO (1 + ^/0^' results in the same frequency shift as ob- 
tained in the previous Appendix, Eq. (IA12I) . 



Appendix C: Qubit probed by tank circuit 

In this Appendix we consider a qubit coupled capacitively 
to the series LCR (tank) circuit [see Fig.Oc)]. The tank cir- 
cuit consists of an inductor Lt and a capacitor Ct, while dis- 
sipation is described by the resistor Rt. The qubit is con- 
sidered to be coupled to the tank circuit through the coupling 
capacitance, which for uniformity we again denote by Cnr 
(even though there is no NR in the scheme considered in this 
Appendix), in parallel to the tank's capacitance Ct- The ef- 
fect of the qubit on the tank circuit can be described by re- 



placing the tank capacitance Ct with Ct = Ct + Cgff , where 
the effective capacitance of the Cooper-pair box is given by 
Eq. (IBll) . The geometric capacitance Cgeom gives only a con- 
stant contribution to the tank capacitance Ct , while the quan- 
tum capacitance Cq <C Cq = Ct + Cgeom is defined by the 
derivative of the average extra Cooper-pair number on the is- 
land (n). 

The tank circuit is biased by the current /b = I a coscjrf^. 
The output voltage is given by Vt = Va cos{uJrft + 6). Then 
from the equation for the voltage we obtain for the phase shift 



tan l9 = Qo 2 



UJo 



Cq 
Co 



(CI) 



1 1 [Lt 

-y==, Acj = cjrf - ojo, Qo = -TT- (C2) 

V^tCo Lit V Co 



The measured value can be either the voltage shift 6 at reso- 
nance frequency (Aco' = 0)^^'^^'^^ 



tan^ = Qo 



Cq 
Co' 



(C3) 



or the resonance frequency shift (at which the voUage shift 



Aw 

too 



Cq 
2Co' 



(C4) 





0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 

FIG. 7: (Color online) LZS interferometry probed via a quantum 
capacitance, (a) The quantum capacitance Cq of the qubit versus the 
energy bias (rig) and the driving amplitude (n^). Arrows show the 
values of and Ug at which the graphs (b) and (c) are plotted as 
functions of rig and n^, respectively. The upper curves were shifted 
for clarity. 
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Both are proportional to the quantum capacitance Cq. 

For the sake of illustration, in addition to Fig. [3l we also 
demonstrate in Fig. |7] the direct LZS interferometry calculated 
for the quantum capacitance for the parameters of Ref. [14]: 
Ejo/h = 12.5 GHz, Ec/h = 24 GHz, cj/27r = 4 GHz, 
ksT/h = 1 GHz, and also we have taken a = 0.005, B = 
0.5. We note that besides the difference in the parameters, 
in Fig. [3] the frequency shift Auj was plotted, while in Fig. [7] 
the quantum capacitance Cq was shown. Both figures were 
calculated by numerically solving the Bloch equation. 

Finally, it is worthwhile emphasizing that for simplicity we 
have assumed that the qubit's dynamics is much faster than the 
resonator's dynamics. In the general case, the cooperative dy- 
namics of the qubit-resonator system should be studied, as e.g. 



in Ref. fs?]. However, a simplification can be made because 
the stationary oscillations in the nonlinear system (either NR 
or tank circuit), influenced by the qubit's dynamics, can be re- 
duced to oscillations in the linear system, as was studied in 
Ref. [17J. In that work, the Krylov-Bogolyubov technique of 
asymptotic expansion was used. This technique describes the 
influence of the qubit as shifts of both the effective damping 
factor and the effective coefficient of elasticity. In analogy to 
the results of Ref. ifTvll . for the system considered here, this 
means that not only the voltage shift is related to the qubit's 
capacitance Cq [see Eq. (ICIH , but also the voltage magnitude 
Va is defined by Cq. This, in particular, explains the experi- 
mental results presented in Fig. 3 by Paila et al. [15.1 . 



* Electronic address: lsshevchenko@ilt.kharkov.ua] 

^ M. R Blencowe, Contemp. Phys. 46, 249 (2005); M. Poot and H. 
S. J. van der Zant, Phys. Rep. 511, 273 (2012). 

^ G. A. Steele, A. K. Huttel, B. Witkamp, M. Poot, H. B. Meer- 
waldt, L. P Kouwenhoven, and H. S. J. van der Zant, Science 325, 
1103 (2009); B. Lassagne, Yu. Tarakanov, J. Kinaret, D. Garcia- 
Sanchez, and A. Bachtold, ibid. 325, 1107 (2009). 

^ M. A. Sillanpaa, J. Sarkar, J. Sulkko, J. Muhonen, and P J. Hako- 
nen, Appl. Phys. Lett. 95, 011909 (2009). 

J. Q. You and E Nori, Nature (London) 474, 589 (2011); Phys. 
Today 58(11), 42 (2005). 

^ E. K. Irish and K. C. Schwab, Phys. Rev B 68, 155311 (2003). 

^ 1. Martin, A. Shnirman, L. Tian, and P Zoller, Phys. Rev B 69, 
125339 (2004); M. R. Geller and A. N. Cleland, Phys. Rev A 71, 
032311 (2005); C. P Sun, L. E Wei, Y X. Liu, and E Nori, Phys. 
Rev A 73, 022318 (2006); S. Ashhab and E Nori, Phys. Rev A 
81, 042311 (2010); Y X. Liu, A. Miranowicz, Y B. Gao, J. Bajer, 
C. P Sun, and E Nori, Phys. Rev A 82, 032101 (2010). 

^ A. D. O'Connell, M. Hofheinz, M. Ansmann, R. C. Bialczak, M. 
Lenander, E. Lucero, M. Neeley, D. Sank, H. Wang, M. Weides, 
J. Wenner, J. M. Martinis, and A. N. Cleland, Nature 464, 697 
(2010). 

^ M. D. LaHaye, J. Suh, P M. Echternach, K. C. Schwab, and M. 

L. Roukes, Nature 459, 960 (2009). 
^ S. N. Shevchenko, S. Ashhab, and E Nori, Phys. Rep. 492, 1 

(2010). 

L. Childress and J. Mclntyre, Phys. Rev. A 82, 033839 (2010); 
L. Chotorlishvili, Z. Toklikishvili, A. Komnik, and J. Berakdar, 
Phys. Rev. B 83, 184405 (2011); S. Miyashita, J. Comp. Theor. 
Nanoscience 8, 919 (2011); J. R. Petta, H. Lu, A. C. Gossard, 
Science 327, 669 (2010); P Plotz and S. Wimberger, Eur. Phys. J. 
D 65, 199 (2011); H. Ribeiro, J. R. Petta, and G. Burkard, Phys. 
Rev. B 82, 115445 (2010); A. Zenesini, D. Ciampini, O. Morsch, 
and E. Arimondo, Phys. Rev. A 82, 065601 (2010); J.-N. Zhang, 
C.-P Sun, S. Yi, and E Nori, Phys. Rev. A 83, 033614 (2011). 
N. B. Delone and V. P Krainov, Atoms in Strong Light Fields, 
Springer Ser. Chem. Phys., Vol. 28, Springer, Berlin-Heidelberg 
(1985). 

M. A. Sillanpaa, T. Lehtinen, A. Paila, Yu. Makhlin, L. Roschier, 

and P J. Hakonen, Phys. Rev. Lett. 95, 206806 (2005). 

T. Duty, G. Johansson, K. Bladh, D. Gunnarsson, C. Wilson, and 

P Delsing, Phys. Rev. Lett. 95, 206807 (2005); E Persson, C. M. 

Wilson, M. Sandberg, and P Delsing, Phys. Rev. B 82, 134533 

(2010). 



M. Sillanpaa, T. Lehtinen, A. Paila, Yu. Makhlin, and P Hakonen, 
Phys. Rev. Lett. 96, 187002 (2006). 

A. Paila, J. Tuorila, M. Sillanpaa, D. Gunnarsson, J. Sarkar, Yu. 
Makhlin, E. Thuneberg, and P Hakonen, Quantum Inf. Process. 
8, 245 (2009). 

V.L Shnyrkov, Th. Wagner, D. Born, S.N. Shevchenko, W Krech, 
A.N. Omelyanchouk, E. Il'ichev, and H.-G. Meyer, Phys. Rev. B 
73, 024506 (2006); V. 1. Shnyrkov, D. Born, A. A. Soroka, and W. 
Krech, ibid. 79, 184522 (2009). 

S.N. Shevchenko, S.H.W van der Ploeg, M. Grajcar, E. Il'ichev, 
A.N. Omelyanchouk, and H.-G. Meyer, Phys. Rev. B 78, 174527 
(2008); S.N. Shevchenko, Eur. Phys. J. B 61, 187 (2008). 
D. A. Garanin and R. Schilling, Europhys. Lett. 59, 7 (2002). 
^'^ M. V. Berry, J. Phys. A: Math. Theor. 42, 365303 (2009). 

D. Burgarth, K. Maruyama, and E Nori, Phys. Rev. A 79, 
020305(R) (2009); S. Ashhab, J. R. Johansson, and E Nori, New 
J. Phys. 8, 103 (2006). 

Y Teranishi and H. Nakamura, Phys. Rev. Lett. 81, 2032 (1998); 
H. Nakamura, Nonadiabatic Transition: Concepts, Basic Theo- 
ries, and Applications, (2nd edn., World Scientific, Singapore, 
2012); B. M. Garraway and N. V. Vitanov, Phys. Rev. A 55, 4418 
(1997). 

E. A. Temchenko, S. N. Shevchenko, and A. N. Omelyanchouk, 
Phys. Rev. B 83, 144507 (2011). 

W D. Oliver, Y Yu, J. C. Lee, K. K. Berggren, L. S. Levitov, and 
T. P Orlando, Science 310, 1653 (2005). 

D. A. Garanin and R. Schilling, Phys. Rev. B 66, 174438 (2002). 
G. Sun, X. Wen, B. Mao, J. Chen, Y Yu, P Wu, and S. Han, Nat. 
Commun. 1, 51 (2010), G. Sun, X. Wen, B. Mao, Y Yu, J. Chen, 
W Xu, L. Kang, P Wu, and S. Han, Phys. Rev. B 83, 180507(R) 
(2011). 

P Huang, J. Zhou, E Eang, X. Kong, X. Xu, C. Ju, and J. Du, 
Phys. Rev. X 1, 011003 (2011). 

T. T. Suzuki and Y Yamauchi, Phys. Rev. A 82, 042709 (2010). 

E. Lucero, M. Hofheinz, M. Ansmann, R. C. Bialczak, N. Katz, 
M. Neeley, A. D. O'Connell, H. Wang, A. N. Cleland, and J. M. 
Martinis, Phys. Rev. Lett. 100, 247001 (2008). 

S. Ashhab, J. R. Johansson, A. M. Zagoskin, and E Nori, Phys. 

Rev. A 75, 063414 (2007). 

L. Du and Y Yu, arXiv: 1012.2917 

G. Johansson, L. Tornberg, V.S. Shumeiko and G. Wendin, J. 
Phys.: Condens. Matter 18, S901 (2006). 
^2 Ya. S. Greenberg and E. Il'ichev, Phys. Rev. B 77, 094513 (2008). 



